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I. INTRODUCTION 



T^ ] The discovery of the accelerated expansion of the present Universe at the end of the 1990s started a new era in 

04 ■ our understanding of the Universe. It became clear that besides visible and dark matter which are gravitationally 

clustered in galaxies and galaxy clusters there should exist, in addition, a fundamentally different form of energy in 

Q ' . the Universe - dark energy, accounting for about three quarters of the energy content of the Universe. This dark 

Q^' energy is not gravitationally clustered, but distributed rather uniformly in space. Since it has a high negative pressure, 

which is rather unusual from the point of view of ordinary matter, it can provide for the current acceleration of the 

{^J)' Universe. 

The key property of dark energy lies in its ability to violate various energy conditions. Of greatest interest for 

-^^ ■ cosmology are the following energy conditions: (i) the strong energy condition, (e + 3p) > 0; (ii) the weak/null energy 

^ I condition (e + p) > (where e and p arc the effective energy density and the pressure of the matter fields filling the 

l/~j . Universe). To provide for the accelerated expansion, it is sufficient to violate the strong energy condition. In turn, 

the failure to satisfy the weak energy condition results in an accelerated expansion of the Universe that is even faster 

^^D , than exponential. 

^^ ' If some form of matter violating these energy conditions docs indeed exist in the present Universe, then it may serve 

(T^ , not only as a mechanism for providing the current acceleration. It may as well be present in compact astrophysical 
^) • objects. Indeed, despite the fact that in describing the evolution of the Universe it is assumed that dark energy is 
Cn ' homogeneously distributed in the Universe, it is possible to imagine a situation where inhomogeneities may arise due 
^7^ . to gravitational instabilities. This could lead to a collapse of lumps of dark energy with the subsequent creation of 
J> ' compact configurations - so-called dark energy stars |ll4la| . 

For dark energy stars, typically at least the strong energy condition is violated. At the present time, investigations 
of dark energy stars are done in two main directions: (i) the construction of stationary configurations to demonstrate 
j^ \ the existence of such objects and their potential stability |l|-KLl|: (ii) the study of the possible formation process of 
such compact objects during gravitational collapse |12h16| . 

Although the true nature of dark energy is currently unknown, several ways have been suggested to model it. The 
simplest variant seems to be the assumption that dark energy is nothing else than Einstein's cosmological constant A, 
which arises as the energy density of the vacuum. However, the well-known "cosmological constant problem" has led 
to the consideration of models where the dark energy evolves in time. Perhaps one of the most developed directions 
here arc theories with various types of scalar fields. 

Scalar fields arise naturally in particle physics including string theory. Using them, a wide variety of scalar-field 
dark energy models have been proposed (for a review, see, e.g., J17l.ll8|). One of the tools developed for modeling dark 
energy arc the so-called phantom/ghost fields. Their important property is that they violate the weak/null energy 
conditions. Based on such fields, in Ref. 119[ a simple cosmological model was suggested, where the phantom dark 
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energy is provided by a ghost scalar field with a negative kinetic term (for the further development of such models, 
see the review (iSj). 

In view of such highly unusual properties of matter for which the weak/null energy conditions are violated, localized 
objects containing such matter should differ strongly from ordinary stars as well. A striking example of the effects of 
the presence of such exotic matter is the possibility to allow for traversable wormholes - compact configurations with a 
nontrivial space-time topology. The term "wormhole" has appeared in the middle of the 1950s in the work of Wheeler 
[20| , who suggested the geometrical model of electric charge in the form of a tunnel connecting two space-time regions 
and filled by an electric field. The Wheeler wormhole was not traversable. But, this idea has stimulated a great deal 
of interest in studying models with a nontrivial space-time topology. One of the most significant achievements in this 
direction is the model of a traversable wormhole suggested in 1988 by Morris and Thorne [2l| as a toy model allowing 
for interstellar travel. Moreover, since recent observational data indicates that matter violating the weak/null energy 
conditions may indeed exist in the present Universe |22| . traversable wormholes have received increasing attention 
ever since as objects which could really exist in nature. 

Providing a nontrivial topology, massless ghost scalar fields were employed in the early pioneering work of Refs. [2J, 
1241 (see also earlier work by Bergmann and Leipnik [25[, where they found a solution for a massless ghost scalar 
field, but not dealing with the question of the wormhole interpretation of the solution obtained) . Further examples 
of configurations with nontrivial topology were found in [2a, [23], based on a ghost scalar field with a Mexican 
hat potential. Interestingly, it was found that this system has only regular, stable solutions for a topologically 
nontrivial (wormhole-like) geometry. In |28l433l | traversable Lorentzian wormholes were further investigated, refining 
the conditions on the type of matter fields that would lead to such space-times. A general overview on the subject of 
Lorentzian wormholes and violations of the various energy conditions can be found in the book of Visser |34l | . 

A recent development concerning such wormhole models is based on the expectation, that such topologically non- 
trivial objects, being traversable, could be bound to ordinary matter satisfying all energy conditions. In |35| . we have 
suggested such a model of a mixed configuration, consisting of a wormhole (supported by a massless ghost scalar 
field) filled by a perfect polytropic fluid. We have shown that there exist static, regular solutions describing mixed 
star-plus- wormhole configurations which possess new physical properties that distinguish them from ordinary stars. A 
preliminary stability analysis performed only for the external region of the configuration suggested that those solutions 
could be stable with respect to linear perturbations. 

The objective of the present paper is to continue the study of the influence of a wormhole on the structure and 
the physical properties of compact stars made otherwise from ordinary matter. To provide a nontrivial topology in 
this model, we here employ a massless ghost scalar fleld as one of the simplest possibilities allowing for a nontrivial 
topology. For the ordinary matter filling the wormhole, we choose neutron matter modeled within the polytropic 
approximation [in the form of Eq. (|13p ]. 

The paper is organized as follows: In Sec. [ill the general set of equations is derived, describing equilibrium configu- 
rations consisting of a massless ghost scalar field and ordinary matter approximated by a polytropic equation of state. 
Here, the boundary conditions are also discussed. In Sec. IIIIl we discuss the physical properties of the solutions, the 
radii, the masses, the binding energies and the pressures. We then present the numerically obtained static solutions 
for two sets of polytropic parameters and consider the issue of their stability within the energy approach. Finally, in 
Sec. llVl we summarize the results obtained and address possible observable effects which may follow from the model 
considered. 

II. DERIVATION OF THE EQUATIONS FOR EQUILIBRIUM CONFIGURATIONS 
A. Lagrangian and general set of equations 

We consider a model of a gravitating massless ghost scalar field in the presence of a perfect fiuid. Our starting 
point is the Lagrangian 
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L = -^^^R--d^^d^'ip + L„,. (1) 



Here, ip is the ghost scalar field, Lm is the Lagrangian of the perfect isotropic fluid (where isotropic means that the 
radial and the tangential pressure of the fiuid agree) which has the form L„i = p J36l |37| . This Lagrangian leads to 
the corresponding energy-momentum tensor 

Tt = {e+ p)u,u'^ - 5'yp - a,^9V + Is^d^^df^^, (2) 



where e and p are the energy density and the pressure of the fluid, u' is the four-velocity. For our purposes, it is 
convenient to choose the static metric in the form 
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(3) 



where v and A are functions of the radial coordinate r only, and dfl^ is the metric on the unit 2-spherc. Also, let us 
present the metric in Schwarzschild-type coordinates which are frequently used in modeling wormholes 
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This parametrization of the metric is employed later on to derive the boundary conditions at the wormhole throat at 
the core of the configuration. 

The (q) and (}) components of the Einstein equations for the metric 1^ and the energy- momentum tensor ^ are 
then given by 
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where a "prime" denotes differentiation with respect to r. 

The field equation for the scalar field is obtained by varying the Lagrangian ([1]) with respect to </?, 
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Using the metric ([3]), this equation is integrated to give 
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where D is an integration constant. 

Not all of the Einstein field equations are independent because of the law of conservation of energy and momentum 
T^.f. ~ 0. Taking the i = 1 component of this equation gives 
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The remaining two Einstein equations are then satisfied as a consequence of the Eqs. ([5]), (|6]), and (J9l). 

Thus, we have five unknown functions: i^. A, (/3, e, and p. Keeping in mind that e and p are related by an equation of 
state, there are only four unknown functions. For these functions, there are four equations: the two Einstein equations 
dn]) and (IHl), the scalar-field equation (|S]), and the equation of hydrostatic equilibrium ©. Using the energy-momentum 
tensor ([2]), the right-hand sides of the Eqs. ([5]) and (|6]) take the form 
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Then, taking into account the expressions for components of the energy-momentum tensor ([2]) 
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and using (|5]), we obtain from © the following equation for hydrostatic equilibrium, 
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(12) 



B. Equation of state 

To model the neutron matter filling the wormhole, it is necessary to choose an appropriate equation of state. The 
choice of the equation of state plays a crucial role in modeling the neutron-star structure. A large number of neutron 
matter equations of state have been suggested, providing different sets of densities and pressures for the neutron stars. 
These include a parametric dependence as in the pioneering work of Oppenheimer and Volkoff 1381 . a polynomial 
dependence as in [39|, an equation of state obtained from field theoretical considerations as in [40[, or a unified 
equation of state (see, e.g., Ref. [4l[ where the analytical representations of unified equations of state of neutron-star 
matter are given). 

Here, we employ a simplified variant for the equation of state, where a more or less realistic neutron matter equation 
of state is approximated in the form of a polytropic equation of state. In particular, we employ the following parametric 
relation between the pressure and the energy density of the fiuid, 

2 , P ; 2 (ch) 
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where ni, is the baryon number density, nj^'^ is some characteristic value of ni,, mi, is the baryon mass, and k and 7 
are parameters whose values depend on the properties of the neutron matter. 

For our purpose it is convenient to rewrite the above equation of state in the form 

&' 1+1/" 2 I /io\ 

p = Kp^ , e^pbc +np, (13) 

with the constant K — kc^{ni^'^ mi,)^^'^ , the polytropic index n = 1/(7 — 1), and pi, = ribmi, denotes the rest-mass 
density of the neutron fluid. 

Setting nib — 1-66 x 10"^"' g and n^^ =0.1 fm~ , we consider below configurations with two sets of values for the 
parameters k and 7: 

(i) k = 0.0195 and 7 ~ 2.34 [431, adjusted to fit the equation of state II of Ref. [4^; we denote this choice by EOSl 
in this paper. 

(ii) k = 0.1 and 7 = 2 |43l|. corresponding to a gas of baryons interacting via a vector meson field, as described 
by Zel'dovich [4J, |4^ (see also Ref. [4a| where relativistic configurations with such an equation of state were 
considered); we denote this choice by E0S2 in this paper. 

Introducing the new variable 9 [43 

Pb = Pbc^ , (14) 

where pbc is the density of the neutron fiuid at the wormhole throat (or the center of the star in the case without a 
wormhole) we may rewrite the pressure and the energy density, Eq. (J13p . in the form 

p = Kpi:'^-e-+\ e = [pb.c^ + nKplf'^e) r . (15) 

Making use of this expression, we obtain for the internal region with ^ Q from Eq. (|12p 

2a(n + l)^ = -[l + a(n+l)0]^, (16) 

dr dr 

where a — KpJJ^ /c^ = Pc/iPbcC^) is a constant, related to the pressure Pc of the fluid at the wormhole throat (or at 
the center of the star in the case without a wormhole) . This equation may be integrated to give in the internal region 
with 6* 7^ the metric function e'^ in terms of 9, 
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(17) 



and e"" is the value of e'^ at the throat where 9 = 1. The integration constant Vc is fixed by requiring that the 
space-time is asymptotically flat, i.e., e'' = 1 at infinity. 



C. Internal set of equations 



We first consider the set of equations in the internal region, where 9^0. Here, the system is characterized by three 
unknown functions: A, 9 and (p. These three functions are determined by the three Eqs. (O, (O, and ([5]), and also 
by the relation (|17p . It is convenient to rewrite these equations by introducing the new function M{r)^ 



^^x _ ^ 2GM{r) 



(18) 

The function M(r) can be interpreted as the total mass within the areal radius r. Thus, for r = r^, where r^ denotes 
the outer boundary of the fluid where 9 = (for more details on this, see Sec. IIIip . we obtain the total mass of 
the configuration within the boundary of the star of radius r^. The limit r 
configuration. 

With this function M{r), Eq. (P yields 



cx) then yields the total mass of the 
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For the spherically symmetric case without a wormhole (corresponding in our case to the absence of the scalar field), we 
need to require the boundary condition M(0) = in order to guarantee regularity at the origin [47|. This corresponds 
to the fact that there is no mass associated with the origin, r = 0. For a wormhole, on the other hand, there exists a 
finite minimal value of the radius, the throat r = tq, and this is associated with a finite value of the mass function at 
the throat, M(ro) ^ 0. 

Now, we introduce dimensionless variables 
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where A has the dimension of inverse length, and rewrite Eqs. (J19p and 



cr(n + l)c2 
in the form 

2 



1/2 




\ + (j[n + \)9 di 

where wc have used expression ([8]) and introduced the dimensionless constant 
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Thus, the internal part (r < r^) of the static configurations under consideration is described by the Eqs. (|2ip and 
22)) together with the scalar-field equation 
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D. External set of equations 



Let us next turn to the external part (r > rt) of the solutions outside the fluid. Asymptotic flatness of the external 
solutions requires that the metric function e'^ tends to one at inflnity. This in turn determines the value of the 
integration constant Vc at the throat. 

To find the external solutions, we start again from the Einstein equations ([5]) and ([6]), and the scalar-field equation 
([5]), taking into account that in the external region there is no ordinary matter, i.e., e = p = 0. This leads to the 
following system of equations. 
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These can be rewritten in terms of the dimcnsionless variables w(^), 4>{£,) and i/(^) as follows, 
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This system still contains the parameter cr as a trace of the influence of the fluid on the external solution due to the 
definition of the dimensionless quantities (pn| . 



E. Expansion at the throat 



Finally, we need to consider the appropriate boundary conditions at the wormhole throat corresponding to the core 
of the configurations 



Let us here briefly comment on the nomenclature used. In contrast to the case of ordinary stars having a center at 
the point r = 0, there exists a finite minimal value of the radial coordinate r = r^ corresponding to the radius of the 
throat in the case of the star-plus-wormhole systems considered here. Moreover, because of the presence of ordinary 
matter, the configurations differ also from simple wormholes. To take this into account, we employ the term "core" 
to describe the "throat" area of such star-plus-wormhole systems. 

In terms of the dimensionless variables introduced above, ^ = ^o denotes the coordinate at the core of the configu- 
ration. Here, the dimensionless density of the fluid is characterized by 



00 = 0(^0) - 1. 



(30) 



This condition corresponds to the fact that at the wormhole throat the density of the fluid is pbc- 

In order to determine the initial value of the function w, we consider that the metric Q satisfles the condition 
(1 — b{r)/r) > throughout the space-time l2l|,[3J]. This implies 



b{ro) = Tq, b'{ro) < 1, and b{r) < r for r > Tq. 



(31) 



Comparing the form of the metric ^ with the form of the metric ([3]) used for the derivation of the Eqs. (PT|) and 
p2|). and taking into account expression (fTS)) . together with the dimensionless variables ((20|) . we find 
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where we have introduced the dimensionless function i?(^) = Ab. From this expression, taking into account the 
relation i?(^o)/Co ~ 1 from (j31|) . we obtain the boundary condition for the function v at the throat 
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This boundary condition implies that, as ^ 
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^0, the following expression vanishes. 
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This in turn means that the coefficient in front of the derivative of the function 6 in Eq. ((22|) goes to zero, which leads 
to a singularity in a generic solution. 

To obtain solutions, which are regular at the throat we consider a Taylor series expansion of the functions v and 9 
in the neighborhood of the point S, — ^o- At a point ^ = ^i close to the throat the expansion for v reads to first order 



v{^i) = wo + 'yi(Ci -Co)- 
Substituting this into Eq. pi|) . we obtain for vi the following expression. 
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Regularity of the solutions of Eq. ((22| is achieved, when we assume that, together with (|34|) . the right-hand side of 
(|22p goes simultaneously to zero. Proceeding from this requirement, we obtain the following expression for D: 
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Substituting this into p6p . we finally obtain 



The Taylor series expansion for the function Q at at the point ^ = ^i reads to first order 
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where we obtain for the coefficient 0i, 
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Nonsymmctric wormholes (with respect to the two asymptotically fiat space-times) could also be obtained, but these 
would satisfy a different set of boundary conditions. 

Thus, a static equilibrium solution is obtained as follows: We start the numerical integration at the point ^ = ^1, 
solving numerically the system of equations (PT|) and ([^ subject to the boundary conditions ([55)) with ([55)1 and ([55)1 
for u, and ([5^ with (|40p for 6*. We then proceed with the integration until we reach the point <^ = ^f,, where the 
function 6 becomes zero [35[. Here, the energy density associated with the fiuid vanishes. The surface bounded by 
the radius ^ = ^b represents the boundary of the neutron matter. Note, that the energy density associated with the 
scalar field is, in general, still finite at this boundary. 

Now, the external part of the solution is sought, starting from the surface of the fiuid at ^ = ^b with the boundary 
conditions vify^ and v(fy^ as determined from the internal part of the solution. Requiring asymptotical flatness of 
the space-time finally allows to determine the value of the integration constant v^ from (jl7|) by requiring e" to be 
equal to unity at infinity. (The values of i^c for the examples shown in Fig. [l]are given in the caption.) Thus, the 
complete solution for the configuration under consideration is derived by matching of the internal fiuid solutions given 
by Eqs. ([2T|) - ([23)) with the external solutions obtained from the system (|27p - (p9|) . 



III. STATIC SOLUTIONS 

In this section, we discuss the numerical solutions of the above sets of internal and external equations and their 
physical properties. 



Radii 



The radial coordinate r describes the areal radius of a sphere with area 47rr^ . The throat radius corresponds to the 
areal radius tq, and the neutron matter of the star is contained within the areal radius r\, (denoted R in the tables). 
The gravitational radius of the system corresponds to the areal radius Vg = 2GM/(? , where M is the total mass. In 
dimensionless coordinates, the areal radius is given in terms of the coordinate ^. Thus, ^0 denotes the throat radius 
and ^b the radius of the neutron fluid. 

Another physically relevant radial coordinate is given by the coordinate ^ associated with the proper radius, which 
gives the distance from the throat. It is defined as follows. 
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Then, the proper radius of the fiuid i?prop is obtained in dimensional variables as i?prop = £,b/A. 



B. Mass contributions 



The system under consideration consists of two parts: the internal region tq < r < rf, and the external region 
Tf, < r < cxj. Correspondingly, the energy density of the system is given by the internal energy density, £i,it obtained 
from the expressions dTU]), ^, (0, (HZ]), and (^0)) 
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and the external energy density, Eoxt, 



text - L^oJext ~ 2 ^4 



The total energy density of the configuration is 

St = ei„te(6 - + i^cxte(e - ih). 

As discussed above, the effective mass M(r) inside a surface with radius r is given by |2ll. [34 
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where the integration constant is determined by 6(ro) = rg, Eq. (JSTJ. Employing dimcnsionless variables ()20p and 
B{£,) [see Eq. ([5^] we obtain the dimcnsionless effective mass M{^), 
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with M{r) = M*M{^), where the quantity 
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fixes the scale of the mass. 

The asymptotic value lim M(r) = M corresponds to the total mass of the configuration, while in dimcnsionless 

r— ^00 

units lim M{^) = A^. For later reference, we now subdivide this expression for the total mass into four dimcnsionless 
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components according to 
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the internal part of the mass of the scalar field 
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and the external part of the mass of the scalar field 
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C. Binding energy 

We start from the total energy E of the system given by 

E = Mc^ = M* f Xth + Xfl + MsRnt + Msfc^t) C^ = f A/th + Affl + Afsfint + Msfcxt ) c^ . (47) 



In order to derive a physically motivated expression for the binding energy of the system, let us consider the relativistic 
continuity equation 

(n,u^).^ = 0, (48) 

with the baryon number density of the neutron fluid ni, and the four- velocity u^. It follows that the neutron particle 
number N is given by 

N = / (7ibM")yq^d3.T = — / pbc^^^r^dr, (49) 

Jro ' "16 Jro 

with M° = e^"'^, and the factor y/—g enters because the natural volume element on the spacelike hypersurfaces is 
needed [4^, \^ . The associated energy of A'' free neutrons is given by 

Efb = Nrubc^. (50) 

For a simple neutron star without a wormhole, the binding energy (B.E.) is defined as the difference of the energy 
of N free particles Efb, Eq. (|5n|). and the total energy E, Eq. (H7|) 

B.E. = Efb - E, (51) 

as discussed, e.g., in Ref. [4^. An analogous definition holds for a boson star composed of N massive bosons. Indeed, 
in order to disperse the particles to infinity, one has to supply precisely this amount of energy to the system. 

For the combined star-plus-wormhole system, on the other hand, we still have to consider how to deal with the 
massless ghost scalar field that supports the wormhole. Thus, let us first address an isolated wormhole made from a 
massless ghost scalar field without any neutron matter. With the above set of boundary conditions, such a wormhole 
has zero energy (respectively, zero mass) [3J]. This agrees with the energy of flat space in the absence of a wormhole. 

Let us therefore consider the binding energy of the star-plus-wormhole system as the difference between the following 
contributions: The energy of TV free particles dispersed to infinity with a (therefore) vanishing ghost scalar field and 
the energy of the combined star-plus-wormhole system. Alternatively, we could consider the difference between the 
energy of TV free particles dispersed to infinity together with the energy of an isolated wormhole and the energy of the 
combined star-plus-wormhole system. Both cases yield the same result, namely the above expression (|5ip . obtained 
without a wormhole. 

D. Pressure 

In the presence of a scalar field, the effective pressure becomes anisotropic (see, e.g., [50| and references therein). 
Indeed, as follows from Eq. ([2]), the radial pressure Pr — ~Tl is distinct from the tangential pressure pt = — Tf = — Tg 
of the system. This allows to introduce the fractional anisotropy /a, 

fa={pr-pt)/Pv (52) 

For the star-plus-wormhole system, the fractional anisotropy is given by 

■^"' -ad''+^ + 1/2 [1 - 2a(n + l)v/^] (j)'^ ' 

Thus, at the core ^o of the configurations the fractional anisotropy is fa{£,Q) — 2 [2a'^{n + 1)^q + l] . Beyond the 
surface of the fluid, i.e., when ^ > £,b, the fractional anisotropy is always equal to two. The above expression also 
implies that the numerator never vanishes for ^ > ^o- The denominator, however, might tend to zero at certain values 
of the radius. This would give rise to a tremendous growth of the fractional anisotropy in the vicinity of these radii. 
A related interesting consequence is the fact that the pressure of the neutron fluid and the total internal pressure 
of the system (which includes the pressure contributions from the fluid and the scalar field) differ. Indeed, for the 
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TABLE I: Characteristics of a set of configurations for EOSl for several values of the dimensionless throat radius ^q. The 
radius of the throat _Rth, the areal radius of the fluid R, the proper radius of the fluid -Rprop and the gravitational radius rg 
(last column) are given in kilometers. The total mass M, the mass at the throat Mth, the mass of the fluid Ma, the internal 
part of the mass of the scalar field Afsflnt, and the external part of the mass of the scalar field Msfext are given in solar mass 
units. 



Pbc, g cm 3 


Rth, km 


R, km 


-f^propi i^rii 


M/Mq 


Mth/Mg 


Mh/Mq 


Msflnt/M0 


Msfext /Mq 


Tg, km 










Without a 


wormhole 










1.0x10^'' 


— 


12.3047 


12.4312 


0.1604 


— 


0.1604 


— 


— 


0.4732 


2.0x10" 


— 


13.3411 


13.6760 


0.4093 


- 


0.4093 


- 


— 


1.2073 


3.0x10" 


— 


13.7197 


14.2885 


0.6680 


- 


0.6680 


— 


— 


1.9704 


4.0x10" 


— 


13.8078 


14.6132 


0.9080 


- 


0.9080 


- 


— 


2.6781 


6.0x10" 


— 


13.5903 


14.8356 


1.2965 


- 


1.2965 


- 


— 


3.8241 


8.0x10" 


— 


13.1647 


14.7824 


1.5647 


— 


1.5647 


— 


— 


4.6152 


1.0x10^^ 


— 


12.6865 


14.6083 


1.7387 


— 


1.7387 


— 


— 


5.1286 


l.SxlO^'^ 


— 


11.5667 


14.0149 


1.9264 


— 


1.9264 


— 


— 


5.6822 


l.SxlO^'^ 


— 


11.0034 


13.6572 


1.9495 


— 


1.9495 


— 


— 


5.7504 


2.0x10^^ 


— 


10.6733 


13.4338 


1.9466 


— 


1.9466 


— 


— 


5.7418 


S.OxlO^'^ 


— 


9.4529 


12.5362 


1.8482 


— 


1.8482 


— 


— 


5.4516 


6.0x10^5 


— 


7.9586 


11.3502 


1.5434 


— 


1.5434 


— 


— 


4.5523 


l.OxlO^'^ 


— 


7.5888 


11.1383 


1.3675 


— 


1.3675 


— 


— 


4.0336 


1.5xlO^« 


— 


7.7045 


11.3870 


1.3082 


- 


1.3082 


— 


— 


3.8588 


3.0x10^® 


— 


8.1252 


12.0257 


1.3375 


- 


1.3375 


- 


— 


3.9452 


e.oxio^'' 


— 


8.2139 


12.2221 


1.3791 


- 


1.3791 


— 


— 


4.0679 


1.0x10^^ 


- 


8.1649 


12.1958 


1.3822 


- 


1.3822 


- 


- 


4.0771 










Co = 


0.1 










l.OxlO" 


0.4306 


12.2483 


12.3651 


0.1567 


0.1460 


0.1562 


-0.1406 


-0.0049 


0.4623 


2.0x10" 


0.4844 


13.2837 


13.6042 


0.3999 


0.1642 


0.3985 


-0.1576 


-0.0052 


1.1795 


3.0x10" 


0.5190 


13.6642 


14.2137 


0.6526 


0.1760 


0.6501 


-0.1680 


-0.0054 


1.9250 


4.0x10" 


0.5450 


13.7564 


14.5375 


0.8871 


0.1848 


0.8831 


-0.1755 


-0.0053 


2.6165 


6.0x10" 


0.5839 


13.5501 


14.7617 


1.2666 


0.1980 


1.2596 


-0.1858 


-0.0051 


3.7361 


8.0x10" 


0.6132 


13.1368 


14.7123 


1.5285 


0.2079 


1.5179 


-0.1925 


-0.0048 


4.5086 


1.0x10^^ 


0.6369 


12.6712 


14.5431 


1.6984 


0.2159 


1.6839 


-0.1970 


-0.0044 


5.0096 


1.5x10^^ 


0.6823 


11.5849 


13.9683 


1.8812 


0.2313 


1.8560 


-0.2026 


-0.0036 


5.5488 


1.8x10^^ 


0.7038 


11.0426 


13.6245 


1.9035 


0.2386 


1.8715 


-0.2035 


-0.0031 


5.6145 


2.0x10^^ 


0.7165 


10.7269 


13.4115 


1.9005 


0.2429 


1.8638 


-0.2033 


-0.0029 


5.6058 


3.0x10^^ 


0.7677 


9.5830 


12.5737 


1.8047 


0.2603 


1.7441 


-0.1976 


-0.0020 


5.3233 


6.0x10^^ 


0.8637 


8.3183 


11.5589 


1.5127 


0.2928 


1.3893 


-0.1684 


-0.0009 


4.4620 


1.0x10^"^ 


0.9420 


8.1305 


11.3823 


1.3321 


0.3194 


1.1630 


-0.1497 


-0.0005 


3.9291 


1.5x10^'' 


1.0092 


8.2145 


11.3384 


1.2135 


0.3422 


1.0299 


-0.1581 


-0.0004 


3.5794 


3.0x10^^ 


1.1354 


7.6987 


10.2505 


0.9529 


0.3849 


0.8780 


-0.3094 


-0.0005 


2.8108 


6.0x10^*^ 


1.2774 


5.6583 


7.4028 


0.7338 


0.4331 


1.7164 


-1.4149 


-0.0009 


2.1643 


l.OxlO" 


1.3933 


4.4538 


5.6845 


0.7137 


0.4724 


5.9667 


-5.7242 


-0.0012 


2.1053 










Co = 


0.2 










l.OxlO" 


0.8612 


12.1022 


12.1922 


0.1492 


0.2920 


0.1472 


-0.2703 


-0.0197 


0.4400 


2.0x10" 


0.9689 


13.1271 


13.4107 


0.3805 


0.3285 


0.3748 


-0.3016 


-0.0213 


1.1222 


3.0x10" 


1.0380 


13.5051 


14.0074 


0.6206 


0.3519 


0.6104 


-0.3200 


-0.0217 


1.8306 


4.0x10" 


1.0900 


13.5984 


14.3220 


0.8431 


0.3695 


0.8275 


-0.3323 


-0.0217 


2.4867 


6.0x10" 


1.1678 


13.3986 


14.5326 


1.2021 


0.3959 


1.1744 


-0.3474 


-0.0208 


3.5457 


8.0x10" 


1.2263 


12.9932 


14.4724 


1.4481 


0.4158 


1.4068 


-0.3549 


-0.0195 


4.2714 


l.OxlO" 


1.2738 


12.5352 


14.2933 


1.6057 


0.4318 


1.5498 


-0.3578 


-0.0182 


4.7363 


1.5x10^^ 


1.3646 


11.4616 


13.6905 


1.7671 


0.4627 


1.6731 


-0.3537 


-0.0150 


5.2122 


1.8x10^^ 


1.4076 


10.9216 


13.3251 


1.7796 


0.4772 


1.6626 


-0.3468 


-0.0134 


5.2492 


2.0x10^^ 


1.4330 


10.6050 


13.0949 


1.7707 


0.4858 


1.6386 


-0.3413 


-0.0125 


5.2228 


3.0x10^^ 


1.5353 


9.4222 


12.1273 


1.6464 


0.5205 


1.4455 


-0.3104 


-0.0093 


4.8561 


6.0x10^^ 


1.7273 


7.6148 


10.1804 


1.2588 


0.5856 


0.9359 


-0.2564 


-0.0063 


3.7130 


1.0x10^*^ 


1.8840 


6.0863 


8.1143 


0.9990 


0.6387 


0.6391 


-0.2722 


-0.0067 


2.9466 


1.5x10^*^ 


2.0185 


4.9563 


6.4447 


0.9069 


0.6843 


0.6076 


-0.3768 


-0.0082 


2.6751 


3.0x10" 


2.2709 


4.0842 


4.7288 


0.9409 


0.7699 


1.6547 


-1.4735 


-0.0101 


2.7754 


6.0x10" 


2.5549 


4.0725 


3.9632 


1.0863 


0.8662 


10.3705 


-10.1402 


-0.0102 


3.2041 


l.OxlO" 


2.7867 


4.3200 


3.6586 


1.2319 


0.9448 


44.7764 


-44.4798 


-0.0095 


3.6336 



TABLE II: Characteristics of a set of configurations for E0S2. The notations are the same as in Table [H 
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Pbc, g cm ^ 


-Rth, km 


R, km 


^prop 1 Kill 


M/Mq 


Mth/Me 


Mh/Mq 


Msflnt/M0 


Msfext/Mo 


Tg, km 










Without a 


wormhole 










l.OxlO^''^ 


— 


35.1484 


35.3940 


0.2777 


— 


0.2777 


— 


— 


0.8192 


3.0x10^^ 


— 


34.0595 


34.7518 


0.7573 


— 


0.7573 


— 


— 


2.2338 


5.0x10^^ 


— 


33.0649 


34.1519 


1.1532 


- 


1.1532 


— 


— 


3.4014 


7.0x10" 


— 


32.1529 


33.5906 


1.4818 


- 


1.4818 


— 


— 


4.3708 


l.OxlO" 


— 


30.9169 


32.8113 


1.8760 


- 


1.8760 


— 


— 


5.5335 


2.0x10" 


— 


27.6628 


30.6535 


2.6482 


— 


2.6482 


— 


— 


7.8112 


3.0x10" 


— 


25.3072 


28.9877 


2.9832 


— 


2.9832 


— 


— 


8.7993 


4.0x10" 


— 


23.5170 


27.6575 


3.1213 


— 


3.1213 


— 


— 


9.2068 


6.0x10" 


— 


20.9771 


25.6722 


3.1564 


— 


3.1564 


— 


— 


9.3101 


8.0x10" 


— 


19.2619 


24.2625 


3.0779 


— 


3.0779 


— 


— 


9.0787 


l.OxlO" 


— 


18.0313 


23.2159 


2.9715 


— 


2.9715 


— 


— 


8.7647 


1.5x10" 


— 


16.1136 


21.5291 


2.7156 


— 


2.7156 


— 


— 


8.0099 


2.0x10" 


— 


15.0592 


20.5768 


2.5160 


— 


2.5160 


— 


— 


7.4212 


3.0x10" 


— 


14.1095 


19.7367 


2.2541 


— 


2.2541 


— 


— 


6.6487 


6.0x10" 


— 


14.0299 


19.8738 


1.9703 


— 


1.9703 


— 


— 


5.8117 


l.OxlO" 


— 


14.8886 


20.9756 


1.9447 


— 


1.9447 


— 


— 


5.7360 


1.5x10" 


— 


15.5199 


21.8035 


2.0028 


— 


2.0028 


— 


— 


5.9076 


3.0x10" 


- 


15.8192 


22.3168 


2.1019 


- 


2.1019 


- 


- 


6.1998 










Co = 


0.1 










l.OxlO" 


1.1374 


35.0540 


35.2774 


0.2717 


0.3856 


0.2708 


-0.3727 


-0.0121 


0.8014 


3.0x10" 


1.1374 


33.9788 


34.6415 


0.7409 


0.3856 


0.7383 


-0.3714 


-0.0116 


2.1855 


5.0x10" 


1.1374 


32.9983 


34.0493 


1.1283 


0.3856 


1.1240 


-0.3701 


-0.0112 


3.3282 


7.0x10" 


1.1374 


32.0990 


33.4947 


1.4500 


0.3856 


1.4441 


-0.3688 


-0.0108 


4.2771 


l.OxlO" 


1.1374 


30.8813 


32.7260 


1.8360 


0.3856 


1.8276 


-0.3669 


-0.0103 


5.4154 


2.0x10" 


1.1374 


27.6798 


30.6016 


2.5927 


0.3856 


2.5764 


-0.3606 


-0.0087 


7.6474 


3.0x10" 


1.1374 


25.3698 


28.9684 


2.9217 


0.3856 


2.8980 


-0.3544 


-0.0075 


8.6180 


4.0x10" 


1.1374 


23.6247 


27.6755 


3.0584 


0.3856 


3.0276 


-0.3482 


-0.0066 


9.0210 


6.0x10" 


1.1374 


21.1631 


25.7582 


3.0959 


0.3856 


3.0518 


-0.3363 


-0.0052 


9.1318 


8.0x10" 


1.1374 


19.5213 


24.4165 


3.0229 


0.3856 


2.9665 


-0.3249 


-0.0043 


8.9165 


l.OxlO" 


1.1374 


18.3622 


23.4391 


2.9230 


0.3856 


2.8551 


-0.3141 


-0.0036 


8.6218 


1.5x10" 


1.1374 


16.6197 


21.9273 


2.6850 


0.3856 


2.5915 


-0.2896 


-0.0025 


7.9197 


2.0x10" 


1.1374 


15.7405 


21.1547 


2.5036 


0.3856 


2.3884 


-0.2686 


-0.0018 


7.3847 


3.0x10" 


1.1374 


15.1024 


20.6263 


2.2768 


0.3856 


2.1276 


-0.2353 


-0.0011 


6.7156 


6.0x10" 


1.1374 


15.7007 


21.3879 


2.0735 


0.3856 


1.8706 


-0.1823 


-0.0005 


6.1160 


l.OxlO" 


1.1374 


17.0156 


22.7755 


2.0881 


0.3856 


1.8611 


-0.1583 


-0.0003 


6.1592 


1.5x10" 


1.1374 


18.1359 


23.8678 


2.1448 


0.3856 


1.9081 


-0.1486 


-0.0002 


6.3264 


3.0x10" 


1.1374 


19.8761 


25.3427 


2.1975 


0.3856 


1.9480 


-0.1359 


-0.0002 


6.4818 










eo = 


0.15 










l.OxlO" 


1.7061 


34.9308 


35.1277 


0.2659 


0.5784 


0.2640 


-0.5492 


-0.0273 


0.7844 


3.0x10" 


1.7061 


33.8677 


34.4974 


0.7252 


0.5784 


0.7193 


-0.5463 


-0.0263 


2.1391 


5.0x10" 


1.7061 


32.8960 


33.9080 


1.1043 


0.5784 


1.0946 


-0.5434 


-0.0253 


3.2572 


7.0x10" 


1.7061 


32.0063 


33.3577 


1.4190 


0.5784 


1.4056 


-0.5406 


-0.0244 


4.1856 


l.OxlO" 


1.7061 


30.8018 


32.5950 


1.7965 


0.5784 


1.7776 


-0.5363 


-0.0232 


5.2990 


2.0x10" 


1.7061 


27.6376 


30.4894 


2.5357 


0.5784 


2.4993 


-0.5223 


-0.0197 


7.4792 


3.0x10" 


1.7061 


25.3584 


28.8739 


2.8559 


0.5784 


2.8033 


-0.5087 


-0.0171 


8.4239 


4.0x10" 


1.7061 


23.6387 


27.5959 


2.9876 


0.5784 


2.9197 


-0.4956 


-0.0150 


8.8123 


6.0x10" 


1.7061 


21.2228 


25.7085 


3.0202 


0.5784 


2.9243 


-0.4706 


-0.0119 


8.9084 


8.0x10" 


1.7061 


19.6208 


24.3938 


2.9446 


0.5784 


2.8236 


-0.4477 


-0.0098 


8.6853 


l.OxlO" 


1.7061 


18.4970 


23.4397 


2.8426 


0.5784 


2.6990 


-0.4266 


-0.0083 


8.3846 


1.5x10" 


1.7061 


16.8313 


21.9741 


2.5989 


0.5784 


2.4079 


-0.3816 


-0.0058 


7.6657 


2.0x10" 


1.7061 


16.0135 


21.2281 


2.4094 


0.5784 


2.1817 


-0.3463 


-0.0044 


7.1067 


3.0x10" 


1.7061 


15.4550 


20.6941 


2.1579 


0.5784 


1.8796 


-0.2972 


-0.0029 


6.3648 


6.0x10" 


1.7061 


16.0124 


21.0850 


1.8223 


0.5784 


1.4805 


-0.2350 


-0.0016 


5.3750 


l.OxlO" 


1.7061 


16.8904 


21.5848 


1.5973 


0.5784 


1.2302 


-0.2101 


-0.0012 


4.7113 


1.5x10" 


1.7061 


16.9757 


21.1650 


1.3620 


0.5784 


0.9774 


-0.1927 


-0.0012 


4.0173 


3.0x10" 


1.7061 


13.5704 


16.5958 


0.8878 


0.5784 


0.4483 


-0.1375 


-0.0014 


2.6186 
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FIG. 1: The metric functions gtt ~ e", Qrr ~ —e , and the total energy density £t from (|42|l and (|43|l (in units of ptcC ) for 
two equations of state EOSl (left panel) and EOS2 (right panel). The graphs are shown for configurations with the maximum 
mass Ad according to Tables HI and HIl The solid and the dashed curves correspond to ^o ~ 0.1 and ^o = 0.2 for EOSl, and 
5o = 0.1 and ^o ~ 0.15 for E0S2, respectively. The thin vertical line marked by ^b corresponds to the boundary of the fluid 
where 6 = 0. The thin vertical lines marked by ^o correspond to the core of the respective configurations. The values of the 
total energy density at the core are for EOSl £t{^o = 0.1) ^ -59.34 and Sti^o = 0.2) « -14.17; for EOS2 et{^o = 0.1) « -82.00 
and et{^o = 0.15) ~ —35.89. To provide asymptotic flatness of the solutions, the value of the constant Uc is chosen for EOSl 
at Co = 0.1 i^c ~ -1.919, at ^o = 0.2 Uc ~ -1.865; for EOS2 at ^o = 0.1 v^ ~ -1.472, at ^o = 0.15 v^ « -1.455. 



neutron fluid the pressure p is given by Eq. (jlSp . In turn, from Eqs. 
radial internal pressure of the system 



Pint 



^-[tU, 



PbcC 



aO"^ 



2~^ 



|), (dni, dUl), and (Uni) in (HH), one finds for the 



l + a{n + l)0 



1 2 



1 + o-(?i + 1) 



(53) 



Taking into account Eqs. (|33p and (|37)) . the ratio of the total radial internal pressure ((53)) and the fluid pressure ([T5 
at the core of the configuration is given by 



Pint 

P 



1 



2a^n + l)Cl 



(54) 



In principle, for large values of ^o f-nd a the absolute value of this ratio can be small. However, for the configurations 
of most interest to us, i.e., for the potentially stable configurations, a is less than one. Then, Eq. (52) shows that 
the absolute value of this ratio is typically considerably greater than one when considering also configurations with 
^0 less than one. Obviously, the difference in these pressures arises because of the presence of the ghost scalar field. 

E. Numerical results 

Examples of numerical solutions of the sets of equations subject to the appropriate sets of boundary conditions are 
presented in Fig. [TJ In particular, we exhibit the metric functions gu = e'^ and grr = — e^ and the distribution of the 
total energy density et of the system consisting of the scalar-field energy density and the fluid energy density. For the 
graphs, we have determined e" from Eq. ([T7)) . and e^ follows from Eq. p^ together with (|20p 



l-2a(7i + l)- 



(55) 
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FIG. 2: The case of EOSl: 

(a) Total mass of the configurations (in solar mass units), M/Mq, versus the core density ptc for various values of the dimen- 
sionless throat radius ^o. Stable configurations should reside to the left of the first mass peak. 

(b) The dimensionless binding energy given by Eq. (|51|) versus the core density phc for various values of the dimensionless throat 
radius ^o- 
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FIG. 3: The case of E0S2: 

(a) Total mass of the configurations (in solar mass units), M/Mq, versus the core density ptc for various values of the dimen- 
sionless throat radius ^o. Stable configurations should reside to the left of the first mass peak. 

(b) The dimensionless binding energy given by Eq. (|51|) versus the core density phc for various values of the dimensionless throat 
radius ^o- 
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(Note, that this expression diverges at ^ = £,0, as discussed in connection with the boundary conditions.) The 
total energy density St is determined according to Eq. (|44p , and consists of the internal energy density (^^ and the 
external energy density (|43l) . As seen from Fig. [TJ at the boundary ^ = ^b of the fluid the total energy density is almost 
vanishing, since the scalar-field "tail" beyond £,b gives only a very small contribution to the total energy density. 

We exhibit a number of characteristic properties of the star-plus- wormhole systems in Tables U and [ll] for the two 
equations of state EOSl and E0S2, respectively. Choosing ^0 = 0, we start in each case with a sequence of ordinary 
neutron stars for later reference. We then consider sequences of star-plus-wormhole systems for several fixed finite 
values of ^o- Each sequence of configurations is obtained by increasing the core density of the fluid pbc in a physically 
relevant range of values. 

The tables then exhibit for each configuration the areal radius of the throat Rth, the areal radius of the fluid R, the 
proper radius of the fluid i?prop and the gravitational radius of the system r^, all in units of kilometers. In addition, 
the total mass M, and the mass contributions Mth, Ma, MgUnt and Mgfoxt are presented in the Tables, employing 
solar mass units. (In order not to encumber the paper, only a few values of S,o have been selected in the Tables.) 

The dependence of the total mass (in solar mass units) on the core density pbc (in grams per cubic centimeter) is 
shown in Figs.[2][a) and[3l^a) for the two equations of state, respectively, and several values of the core radius ^o- The 
corresponding values of the binding energy Eq. (j5ip are shown in Figs. [DJb) and|31Jb). In these figures, the curves 
labelled by "wo WH" correspond to the configurations without a wormhole (and correspondingly without a scalar 
field), representing ordinary neutron stars modeled by the equation of state (J13p . 

We note that as a function of the core density p^c, the total mass M of the configurations rises monotonically to 
a maximum and then decreases again. This is typical for this type of configurations. In Tooper's paper [46[, for 
instance, ordinary neutron stars were investigated in detail for various values of the polytropic index n. In particular, 
it was shown that the first peak in the mass corresponds to the point dividing stable and unstable neutron-star 
configurations. This first mass peak is reached at a critical value of the core density, p^J ■ 

For the equations of state EOSl and E0S2 used here, the values of n are n ~ 0.75 and n = 1, respectively. Here, 
the typical behavior of the mass function of the case without a wormhole persists in the presence of the wormhole. 
Thus, while the first mass peak persists, the inclusion of a wormhole leads to a decrease of the height of the mass 
peak with increasing wormhole throat. This decrease is stronger for EOSl than for E0S2. 

A necessary, but not sufficient, condition for the stability of the configurations is the positivity of the binding 
energy. Configurations with a negative binding energy are unstable against dispersal of the matter to infinity. In the 
case without a wormhole considered in [46|, the binding energy always has a positive first peak for n <3. Moreover, 
it remains always positive to the left of the peak, and becomes negative to the right of the peak at some large value of 
Pbc- Configurations to the left of the first peak are stable, while configurations to the right are unstable, because it is 
energetically favorable for them to make a transition to a state with the same particle number N, but with a smaller 
central density pbc- 

For the star-plus-wormhole systems, we expect an analogous behaviour concerning their stability. To the left of 
the peak, the solutions should be energetically stable, while to the right of the peak instability would set in. This 
conclusion is supported by a stability analysis via catastrophe theory, as often applied to boson stars (see, e.g., |5ll - 
|53|). Here, one selects an appropriate set of behavior variahle(s) and control parameter(s). In our case, we could 
choose the size of the wormhole throat ^o and the mass M as two control parameters, and the core density pbc as 
a behavior variable. According to catastrophe theory, the stability with respect to local perturbations then changes 
only at turning points, where, in this case, dM jdpbc = for fixed ^o- Thus, stability should change at the maximum 
of the mass. However, a conclusive answer on the question of stability should be obtained from a detailed stability 
mode analysis of the star-plus-wormhole systems. 

IV. CONCLUSIONS 

In this work, we have continued the investigations of star-plus-wormhole configurations suggested in [35|- These 
objects have a nontrivial wormhole-like topology with a tunnel filled by ordinary matter. For obtaining the nontrivial 
topology, we have used one of the simplest forms of matter violating the weak/null energy conditions - a massless 
ghost scalar field. The matter filling the wormhole has been taken in the form of neutron matter approximated by an 
equation of state of the form P^ . 

Since the resulting star-plus-wormhole configurations possess properties both of wormholes and of ordinary stars, 
our goal has been to clarify the influence of the presence of a ghost scalar field on the physical characteristics of such 
well-studied configurations as neutron stars. Our main emphasis has been (i) the study of the dependence of the mass 
of the mixed star-plus-wormhole configurations on the core density and on the size of the wormhole throat; (ii) the 
study of the binding energy and the pressure of the mixed star-plus-wormhole configurations. 

Our results can be summarized as follows: 
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(1) For two polytropic equations of state denoted by EOSl and E0S2 [see Eq. ([T3l) below], sequences of static, 
regular solutions describing neutron-star-plus-wormhole systems have been constructed numerically, by solving 
the coupled Einstein-matter equations subject to appropriate boundary conditions. Examples of such solutions 
are shown in Fig. [TJ 

(2) For the above equations of state, the physical properties of the configurations are characterized by two parame- 
ters - the core density of the fluid and the dimensionlcss throat radius. We have evaluated the total mass of the 
respective configurations, and also the mass contributions of the individual components the system consists of. 
Examples of their values are given in Tables HI and HIl Likewise, the tables exhibit the values of various physically 
relevant radii of these configurations. 

(3) As a function of the core density of the fluid, the total mass of the star-plus-wormhole configurations rises 
monotonically to a maximum and then decreases again, as is typical also for ordinary neutron stars. Since 
neutron stars located to the left of the maximum are stable, we have conjectured that the same may hold for 
the star-plus-wormhole configurations constructed here. 

(4) In the region where stability may be possible, the masses of configurations with a wormhole are always less than 
the masses of ordinary neutron stars without a scalar field. Correspondingly, the maximum values for the mass 
of star-plus-wormhole configurations are always less than for configurations without a wormhole present, as seen 
in Figs. Hand El 

The differences found between the characteristic properties of neutron matter configurations with and without a 
ghost scalar field, of course, are not the only ones that are possible. Further examples may include the following 
phenomena: 

(i) The presence of an intrauniverse wormhole may lead to the fact that a distant observer will see two stars 
associated with the two mouths of the wormhole, having the same or very similar characteristics. 

(ii) Since the neutron matter can, in principle, move freely through the tunnel, this may lead to oscillations of 
matter relative to the core of the configuration. As a result, at various moments of time there will be different 
quantities of matter near the two mouths of the wormhole. From the point of view of a distant observer, this 
will look like periodical changes in the mass of the stars and fluctuations of their luminosity. 

(iii) The neutron-star-plus-wormhole configurations considered here might possibly possess favorable conditions to 
convert neutron matter into quark (or strange) matter, conjectured to be present in quark stars [S^ - ISTI . The 
quark matter would not be able to decay into another more energetically favorable state, since it is assumed to 
possess the highest binding energy possible. In this case, it might be possible that the final stage of evolution 
of such configurations will be quark star-plus-wormhole configurations. Such systems should have a number of 
distinctions as compared to neutron-star-plus-wormhole configurations |57| . 

(iv) A further possible observationally interesting effect of the presence of dark energy in mixed neutron-dark-energy 
configurations consists in possible changes in the gravitational wave spectrum as considered in Ref. |ll| . The 
investigations performed there for configurations with a trivial topology can be also done for systems with a 
nontrivial topology as the ones studied here. The consideration of these and similar issues will be the object of 
further studies. 

Let us finally come back to the question of stability of the st ar-p lus-wormhole systems. It is well known, that an 
isolated wormhole with a massless ghost scalar field is unstable (58l - l61j . When filling such a wormhole with neutron 
matter, a stable configuration may possibly result. However, if a careful mode analysis will reveal that the instability 
persists, one may turn to other, stable types of wormholes when considering star-plus-wormhole systems. Possible 
candidates for such wormholes are, for instance, the recently constructed dilatonic Einstein-Gauss-Bonnet wormholes 

MM- 
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